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© Dispersion-compensating diffraction grating. 



© A transmission grating is used to reduce chro- 
matic aberration in an optical system, wherein the 
system's chromatic dispersion is offset and substan- 
tially canceled by the grating's diffraction-induced 
dispersion. The grating comprises a Fresnel-type 
pattern of microscopic facets molded directly into 
one of the system's refractive optical surfaces. The 
facets would typically have a profile height of around 
4*10"~ 5 inch and a profile width of at least 10~ 3 inch. 
In its primary intended application, the invention 
would function to improve the optical performance of 
a Fresnel lens used to concentrate direct sunlight 
into a piped daylight illumination system. For this 
application, about a ten-fold reduction in chromtic 
^dispersion would be attained with only about a 2% 
loss in optical efficiency. 
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© Dispersion-compensating diffraction grating. 

© A transmission grating is used to reduce chromatic aberration in an optical system, wherein the system's 
chromatic dispersion is offset and substantially canceled by the grating's diffraction-induced dispersion. The 
grating comprises a Fresnel-type pattern of microscopic facets molded directly into one of the system's 
refractive optical surfaces. The facets would typically have a profile height of around 4M0" 5 inch and a profile 
width of at least 10~ 3 inch. In its primary intended application, the invention would function to improve the optical 
performance of a Fresnei lens used to concentrate direct sunlight into a piped daylight illumination system. For 
this application, about a ten-fold reduction in chromtic dispersion would be attained with only about a 2% loss in 
optical efficiency. 
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BACKGROUND OF THE INVENTION 



Field of the invention 



This invention pertains to the use of a diffraction grating for correcting chromatic dispersion in an optical 
system. The primary intended application for the invention would be to improve the optical performance of 
a Fresnel lens employed as a concentrator element in a high-concentration solar collector, which would be 
used in a piped daylight illumination system. Other potential applications may include optical instruments 
25 and imaging systems that operate over a limited spectral band. 



Description of the prior art 



30 



Chromatic dispersion in an optical system is generally minimized by means of an achromatic doublet 
comprising two cemented or closely-spaced lens elements composed of optical materials with different 
dispersive powers. By appropriately balancing the doublet's opticaJ power between the two elements, the 
combination can be made precisely achromatic at two specific design wavelengths. Color correction at 
35 more than two wavelengths can be achieved by using three or more glass materials in the optical system. 
Such systems exhibit excellent achromatic imaging performance/but this method of dispersion compensa- 
tion is not practical for many applications due to the high material cost and weight associated with the lens 
elements and complications associated with multiple-element lens systems. In particular, the method is 
unsuitable for Fresnel lens systems such as concentrating solar collectors. 
40 An optical system can also be achromatized by means of a dispersion-compensating diffraction grating 
whose diffraction-induced chromatic dispersion offsets and substantially cancels dispersion induced by 
other elements in the system. Two rudimentary systems of this type are known in the prior art: Foster 
(1958) describes a mechanism comprising two diffraction gratings that function conjunctively to minimize a 
prism's refraction-induced dispersion. Lohmann (1970) describes the use of a dispersion-compensating 
45 grating to offset diffraction-induced dispersion produced by a phase-modulating grating in a phase contrast 
microscope. Dispersion-compensating mechanisms such as these are only applicable to optical configura- 
tions in which the dispersion-producing and dispersion-compensating elements both function merely to 
deflect a coliimated incident beam into a transmitted beam that is also collimated. More general types of 
dispersion, such as that produced by refraction at a lens surface, cannot be corrected by these means. 
50 Furthermore, these systems have the additional disadvantage that the dispersion-compensating mechanism 
requires extra optical elements which add to the system cost and complexity and reduce its optical 
transmittance (due to multiple surface reflection losses). 



SUMMARY OF THE INVENTION 
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An objective of my invention is to provide a means for offsetting chromatic aberration in a general 
optica! system, wherein no specific restriction is placed on the type or number of optical elements whose 
chromatic dispersion can be corrected by a single dispersion-compensating element, and wherein the 
beam's wavefront shape in the vicinity of the dispersion-producing and dispersion-correcting elements need 
not be restricted to the planar form of a collimated beam. 

Another objective of the invention is to achromatize an optical system by a means that does not require 
any additional optical components or surfaces in the system. 

A third objective of the invention is to achieve achromatic imaging in an optical system while 
maintaining the highest possible optical efficiency. 

These objectives are achieved by means of a dispersion-compensating element comprising a blazed 
transmission grating formed as a relief pattern on the surface of one of the system's refractive optical 
elements (e.g.. a lens or prism). In addition to minimizing chromatic dispersion, the grating also contributes 
a marginal amount of optical power to the refractive element. The grating's design configuration is 
determined so that the element simultaneously generates an output beam of a specified wavefront shape 
and compensates for chromatic dispersion at all optical interfaces intercepted by the beam. The grating's 
dual role as a beam-shaping and dispersion-compensating element obviates the need for any additional 
dispersion-compensation components, since, even without dispersion compensation, a refractive surface 
would need to be used in place of the grating to provide the beam-shaping capability. 

FIG. 1 illustrates a typical dispersion-compensated lens element's profile geometry (i.e., its cross- 
sectional geometry in a plane transverse to the grating lines). The grating element comprises a surface 
relief pattern molded directly into the lens surface 1. Its geometric form is similar to that of a curved Fresnel 
lens, but with microscopic facet dimensions. (The magnified portion of FIG. 1 illustrates the individual facet 
elements 2.) The grating lines on an axially symmetric lens would comprise concentric circles with the line 
density varying radially in approximately direct proportion to distance from the lens axis. Near the edge of a 
dispersion-compensated lens with a focal ratio of one, the line density would typically be around 1000 lines 
per inch, which is over an order of magnitude lower than line densities required for most conventional 
grating applications. (The line density is so low because only a small amount of diffractive dispersion is 
required to balance the lens's intrinsic refractive dispersion.) The linear-ramp, Fresnel-type facet structure 
maximizes diffraction efficiency in the primary diffracted order, which would be a 1st or -1st order. (A 
grating with this type of facet geometry is termed a "blazed" grating.) The facet profile height would 
typically be around 4°10 -5 inch over the full lens aperture. (For a line density of 1000 lines per inch, the 
facet height-to-width ratio would be around 1/25.) 

The primary application for which dispersion-compensating gratings are currently being considered 
would be to improve the optical performance of Fresnel lenses used as solar energy concentrators. FIG. 2 
illustrates the profile geometry of a dispersion-compensated Fresnel lens. The lens surface 3 has a two- 
tiered facet structure comprising large-scale (e.g., 25 line-per-inch) Fresnel facets 4 on which small-scale 
(e.g., 1000 line-per-inch) grating facets 5 would be superimposed. The Fresnel facets would provide most of 
the lens's optical power, and the grating structure would provide some additional optical power, in addition 
to compensating for chromatic dispersion. 

Chromatic dispersion between two selected design wavelengths can be eliminated by means of a 
grating structure formed on a single lens surface. The residual dispersion over the visible spectrum would 
typically be lower than that of an uncorrected lens by around a factor of 10. For solar energy applications, 
this level of dispersion correction would make it possible to achieve a geometric concentration of up to 
around 10.000X in a Fresnel lens's focal plane. Diffraction efficiency would be quite high if the system uses 
only a limited portion of the solar spectrum. Under their currently-envisaged application as solar concentra- 
tors in piped illumination systems, dispersion-compensated Fresnel lenses would exhibit spectrally-averag- 
ed diffraction efficiency losses of only about 2% in the primary order, based on a typical air-mass-one solar 
illuminance spectrum. For concentrator photovoltaic applications the efficiency loss would be greater, 
perhaps as high as 10%, depending on the characteristics of the photovoltaic cell; and for solar thermal 
applications that utilize the full solar spectrum the efficiency loss would be around 30%, based on a typical 
air-mass-one solar irradiance spectrum. (The diffraction efficiency loss represents energy that is diffracted 
into extraneous diffracted orders. This energy would form a halo surrounding the sun image in the lens's 
focal plane and would be available for capture, but at a comparatively low concentration level.) 

Since a dispersion-compensating grating would have an extremely low line density and a very low- 
profile facet geometry in relation to conventional diffraction gratings, the structure could be produced by 
standard lens molding processes. For a Fresnel lens, no modification of the fabrication process would be 
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required - the grating structure would simply be incorporated in the molding tool design. The molding tool 
would typically be formed by a precision turning operation, using a straight-edge, single-crystal diamond 
cutter to define the grating facet surfaces. 

Several features distinguish the invention from the previously-cited prior art. The grating element in my 
invention is formed as an integral part of a refractive optical element that functions as both a beam-shaping 
and dispersion-compensating element, whereas in the inventions of Foster (1958) and Lohmann (1970) the 
grating forms a separate and distinct optical element which has no refractive optical power. The blazed 
grating structure used in my invention avoids the 1.V intensity loss of Lohmann's system. Foster's 
invention specifically requires two diffraction gratings, and does not anticipate a dispersion-compensating 
mechanism comprising a single grating element. 

Although my invention requires only one grating surface, this advantage would not preclude its use in 
combination with one or more other such elements in a series arrangement. An essential difference 
between a two-grating system of this type and Foster's two-grating system is that, whereas Foster's 
invention requires two gratings to achieve two-wavelength achromatization. the extra grating in a combina- 
tion system employing my invention would merely function to minimize residual dispersion* (e.g., two 
grating surfaces would have sufficient degrees of freedom to achieve four-wavelength achromatization of a 
specific incident wavefront) or to migrate other optical aberrations. Only one grating would be required to 
achieve two- wavelength achromatization. Combination systems employing more than one dispersion- 
compensating grating will later be given some brief attention, but a general treatment of such systems 
would be beyond the scope of this specification and combination systems are not claimed under this patent 
(although the individual dispersion-compensating grating elements in such a combination are claimed). 

C Note: Foster (1958) defines "residual dispersion" in a manner that is inconsistent with its use here. In 
accordance with convention, "residual dispersion" is defined here as the remaining dispersion in an 
achromatic system's output beam, whereas Foster defines it as the dispersion in the portion of the beam 
between the dispersion-producing and dispersion-compensating elements.) 

The dispersion-compensating means discussed in this specification can be extended in an obvious 
manner to reflection gratings as well as transmission gratings, but only the latter type are claimed under this 
patent. 



BRIEF DESCRIPTION OF THE DRAWINGS 



FIG. 1 illustrates a typical dispersion-compensated lens's profile geometry (i.e.. its sectional geome- 
try in a plane transverse to the grating lines). 

FIG. 2 illustrates the profile geometry of a dispersion-compensated Fresnel lens. 

FIG. 3 is a sectional view of a uniform grating structure under illumination by a plane wavefront. 

FIG. 4 illustrates a derivation of the grating equation based on the phase-matching principle. 

FIG. 5 illustrates a derivation of the bla2e condition based on the phase-matching principle. 

FIG. 6 defines the coordinate system used in deriving the grating efficiency formula. 

FIG. 7 illustrates an application of the phase-matching principle to the simple problem of designing a 
conventional refracting lens surface. 

FIG. 8 illustrates the structural elements of a dispersion-compensating diffraction grating. 

FIG. 9 illustrates a procedure by which a diffracted beam's phase distribution can be determined 
from a specification of its phase over the grating substrate. 

FIG. 10 illustrates a structural distinction between a dispersion-compensating element and the 
refractive optical element on which it is formed. 

DETAILED DESCRIPTION OF THE INVENTION 



Uniform grating configurations 



Prior to developing the theory of dispersion-compensating gratings, we will first review some basic 
principles relating to diffraction gratings in general. FIG. 3 illustrates the configuration that we will initially be 
considering. The grating comprises a set of parallel grooves (or "lines"; e.g., groove 6, shown cross- 
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sectionally in FIG. 3) on a surface separating two optical media whose refractive indexes n are indicated as 
n 0 and n 0 '. We refer to an imaginary surface 7 tangent to the grooves as the "grating substrate". We will 
initially restrict our attention to a "uniform" grating structure, which is characterized by the following 
properties: The substrate is flat; the groove structure has translational symmetry in a direction parallel to the 
5 grating lines; and the structure has periodic symmetry in the direction perpendicular to the grating lines and 
parallel to the substrate. In addition, we shall also assume that the incident beam is collimated (i.e., its 
wavefronts are planar). 

Although a dispersion-compensating grating would typically be nonuniform, a uniform grating model 
would generally fairly accurately represent the diffraction geometry over any small-scale region on the 

io grating aperture over which the grating structure is substantially uniform (i.e.. a region that is large in 
comparison to the grating period, but small in comparison to the grating aperture dimension and the 
curvature radii of the substrate, the grating lines, and the incident wavefront). 

We will initially only consider diffraction geometries in which the incident beam's propagation direction 
is perpendicular to the grating lines. (The analysis will later be extended to include the more general 

75 "conical diffraction" case.) FIG. 3 illustrates an incident plane wavefront 8 impinging on the grating in the 
index-rio medium at an incident angle of 0 relative to the grating substrate. The grating will scatter the beam 
into a number of collimated diffracted beams, one of which is illustrated in FIG. 3 as a plane wave 9 
transmitting into the index-n 0 ' medium with an exit angle e relative to the substrate. 

The diffracted orders' exit angles d (FIG. 3) can be determined by a phase-matching condition: 

20 Considering two reference wavefronts, one in the incident beam and one in a diffracted beam, the optical 
path length between these two wavefronts along any two rays traversing corresponding points on adjacent 
grooves (e.g., adjacent troughs) must differ by some integer number of wavelengths. This is illustrated in 
FIG. 4: Ray segments 10 and 11 traverse the grating from wavefront 12 to wavefront 13. crossing adjacent 
trough points 14 and 15. The portion of segment 11 in the index-n 0 medium is longer than that of segment 

25 10 by a difference L = A sinfl (where A is the grating period) and the portion of segment 10 in the index- 
no' medium is longer than that of segment 11 by a difference l' = a sine'. Thus, the phase-matching 
principle implies that 

30 A sine' . A sinO + 

for some integer m, where X D and \ Q are the wavelengths in the index-n Q and index-n D ' media, respectively. 
X 0 and \ 0 ' are each smaller than the free-space wavelength X by a factor of the refractive index: 
X 0 = X/n 0 ,X o ' = X/n 0 '. (2) 

Equation (1) is the "grating equation" for the m-th diffracted order. In general, a separate diffracted 
order will be generated in the index-n Q ' medium for each integer m. (Corresponding reflected orders will 
also be scattered back into the index-n 0 medium, but these are not of interest for the present application.) 
We will only be concerned with first-order diffraction since a dispersion-compensated optical system would 
generally utilize only the first order, so we will set m equal to 1 in equation (1): 
m = 1. (3) 

With the above substitutions, the grating equation (1) for the first order reduces to 
n 0 ' sinfl' = n 0 sinfl + X/A (1st order). (4) 

In some instances a dispersion-compensating grating's primary diffracted order will actually be a -1st 
order, but in this situation we can alternatively represent the primary order as a + 1st order by defining A as 
a negative quantity. Referring to equation (1). a grating operating in the -1st order (m = -1) with A > O is 
formally equivalent to one operating in the + 1st order (m = +1) with A < O; hence we can set m to + 1 in 
either case and can use the sign of A to indicate whether the grating is actually operating in the + 1st or 
-1st order. 

According to the phase-matching condition described above, the radiance emitted from corresponding 
points of the grating profile on different grooves will converge in phase at an image point where the m-th 
order is brought to a focus. The radiance originating from different parts of the groove profile may combine 
out of phase, however, resulting in low diffraction efficiency in the m-th order. In order to achieve high 
efficiency in a particular diffracted order, the grating profile shape would also be defined by a phase- 
matching condition: Again considering two reference wavefronts, one in the incident beam and one in a 
diffracted beam, the optical path length between these two wavefronts along any two rays traversing the 
groove profile should be identical. As illustrated in FIG. 5, this condition will hold if the grating facets have a 
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lmear-ramp profile shape with an appropriately chosen facet angle a. Considering two ray segments 16 and 
17 that cross facet points 18 and 19 while traversing from wavefront 20 to wavefront 21, the portion of 
segment 17 in the index-n 0 medium is longer than that of segment 16 by a difference L = D sin<e + a ), 
where D is the distance between points 18 and 19. and the portion of segment 16 in the index-n G ' medium 
5 is longer than that of segment 17 by a difference L = D sin(*' + «). Diffraction efficiency will be maximized 
when the phase-matching condition. L \ 0 = L '\ 0 ' holds, or equivalently, 
n 0 sin(<9 +a) = n o sin(0 + <*). (5) 

(Note that this relation is simply Snell's Law, applied at the facet surface.) 

A grating with a linear-ramp facet profile geometry of the type illustrated in FIG. 5 is termed a "blazed 
70 transmission grating"; equation (5) is termed the "blaze condition"; and the particular wavelength X that 
satisfies equation (5) (with $ implicitly defined by the 1st-order grating equation (4)) is termed the "blaze 
wavelength". The facet angle a that satisfies the blaze condition for a particular diffraction geometry and 
blaze wavelength can be obtained from equation (5), which can be equivalently stated 



75 tana * 



20 



A(n 0 cos6 - n Q ' coa6M ' (6) 



Dispersion compensation 



Given a specification of an incident beam's angle of incidence d, its wavelength X. and the exit angle B 
25 defining the direction into which the beam is to be diffracted, a grating structure that will diffract the beam 
efficiently into a first diffracted order propagating in the specified direction can be determined a two-step 
process: We first determine the grating period a so that the diffracted beam has the specified propagation 
angle 6 : and we then determine the grating's facet angle a so that the diffraction efficiency in the first order 
is maximized. (A is determined by the grating equation (4) and a is defined by the blaze condition (6).) 
^ The above design procedure presumes that the grating substrate's orientation relative to the incident 

beam is specified. If the substrate orientation is not predetermined, however, we have an extra degree of 
freedom which can be taken advantage of so that two specified wavelengths can be made to simultaneously 
diffract in the designated direction (taking into consideration the dependence of n Q and n Q ' on wavelength). 
We will first consider how this can be accomplished for two very close wavelengths X g and X g + 5X, where 
^ S\ is a small quantity. (If we interpret <5X as an infintesimal, the exit angle's derivative with respect to 
wavelength. 69 /dX, will be zero at the design wavelength Xg.) In order to obtain good dispersion- 
compensating performance over the visible spectrum the design wavelength Xg would typically be chosen 
in the green part of the spectrum (hence the subscript "g"). 

Before developing the grating design, we will make a couple of generalizations: We suppose that the 
two design wavelengths originate from the same source, but that prior to intercepting the grating they may 
have transmitted through preceding dispersive elements that may have caused their incident angles 9 at the 
grating substrate to diverge slightly. (If the grating is formed on a lens element's back surface, for example, 
refraction at the front surface would induce chromatic dispersion in the beam before it intercepts the 
grating.) The incident angles associated with wavelengths X g and X g + 5X will be indicated as 9 and 9 + 89, 
respectively. We will also suppose that the two design wavelengths are to be brought to a common focus! 
but that after traversing the grating they may need to diverge slightly in order to offset the effect of 
chromatic dispersion at optical interfaces following the grating. (If the grating is located on a lens element's 
front surface, for example, the grating would need to induce some chromatic dispersion in the beam in 
order to offset refraction-induced dispersion at the back surface.) The exit angles for wavelengths X g and X 
+ 5X will be indicated as 6 and 6 + 89 . respectively. 9 9 

The incident and exit dispersion angles 6$ and 89 are specified, but since the grating's substrate 
orientation is not predetermined 9 and 9 are not initially specified. Instead, we specify the deflection angle 

* between the incident and exit beams at wavelength X g . * is the difference between e and 9- 

* = B - 6. (7) 

55 , From the gating equation (4) the diffraction geometry for wavelength X g is constrained by the relation 
n 0 sme = n Q sine + x^ A (8) 

where n 0 and n G are the refractive indexes at X = X g on the grating's incident and exit sides, respectively. 
Indicating the corresponding index values at X = X g + 8\ as n Q + 8n Q and n G ' + Sn 0 \ we similarly obtain 
(n Q + 5n G ) sin (9 + 89 ) = (n G + 5n Q ) sin (6 + 89) + (Xg + SX)/A. (9) 
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Subtracting (8) from (9) and using first-order differential approximations, we obtain 
(5n 0 ' sine' + n Q ' 66 cose') = (<5n Q sine + n 0 SO cose) + 5X,A. (10) 

Using equations (7), (8), and (10) we can solve for the three unknowns 6, e\ and A. Using (7) to 
5 eliminate 9 in (8) and (10) and then eliminating A between (8) and (10) and solving for 0, we obtain 
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(The ratios of infintesimal quantities in (11) may be interpreted as derivatives; e.g.. Sn 0 S\ = dn^dXjX = \ g .) 
Equation (11) determines the grating substrate's orientation relative to the incident beam; and given e we 
then obtain Q and A from (7) and (8). (This procedure may yield a negative value for A, indicating that the 
75 grating is actually operating in the -1st order.) In order to maximize diffraction efficiency for wavelengths 
near Xg, we may define the grating facet angle a by equation (6), with all parameters evaluated at the 
design wavelength Xg. 

The above design procedure results in excellent dispersion correction for narrow-band applications; but 
when designing for a wide spectral band better performance can be obtained by using two widely separated 

20 design wavelengths X r and X b (which would typically be in the red and blue parts of the spectrum, 
respectively). We will indicate the incident angles for wavelengths \ r and \ b as e r and e bi respectively; and 
the corresponding exit angles will similarly be indicated as a/ and e b \ We specify the following design 
parameters: the dispersion angle Ae between the red and blue incident rays, 
Ae = e r - e b , (12) 

25 and the red and blue rays' respective deflection angles * r and * b : 

*r = Qr ; B r , (13) 

*b = *b - *b. (14) 

Applying the grating equation (4), we obtain 

n r sine, = n r sine r + \ r L, (15) 

30 n b ' sine b ' = n b sine b + X^'L; (16) g 

where n r and n/ are the refractive index values (previously indicated as n Q and n Q ) at wavelength x r on the 
grating's incident and exit sides, respectively, and n b and n b ' are the corresponding index values for 
wavelenngth X b . The grating configuration is determined by solving the five equations (12)-(16) for the five 
unknowns e r , e b , e r \ 0 b \ and A. Using (12)-(14) to eliminate e r , 0 r \ and e b ' in (15) and (16) and then 

35 eliminating A between (15) and (16) and solving for e bl we obtain 

tan^b = 

x r (nb' 3in<^fa>) - Xfa <n r " sin<A& * - n r 33.n<A3)) 

40 Cn b " n b' coa< > * Xb <nr' coa<A8 * ^ r ) - n r co©< a$> > . 

<17> 

Given e b , we can obtain e b ' from (14) and can then determine A from (16) ; and 0 f and 0,' are determined 
from (12) and (13). 

When designing a grating for wide-band applications, the facet profile shape would be designed to 
maximize some measure of the grating's spectrally-averaged diffraction efficiency (using an appropriate 
spectral weighting function, such as the solar illuminance spectrum for daylighting applications). The 
straight-profile blazed grating profile shape described above would only be strictly optimum when the 
grating is designed to maximize efficiency at a single design wavelength; but for wide-band applications the 
optimum degree of profile curvature and the performance advantage gained by using a curved facet profile 
would be insignificant, so the grating could in practice have a simple straight profile structure with its blaze 
wavelength chosen to maximize the band-averaged efficiency. (A blazed grating's diffraction efficiency as a 
function of incident angle and wavelength can be predicted using a general efficiency formula, equation 
(35), which will be developed below). Once the grating's substrate orientation and period are determined 
(from equations (17), (14) and (16)), the facet angle a may be defined by choosing an appropriate blaze 
wavelength near the middle of the spectral band, determining this wavelength's exit angle e' from the 
grating equation (4), and then using the blaze condition (6) to calculate a. 
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Diffraction efficiency 

In order to optimize a dispersion-compensating grating's design and fully characterize its performance 
we must be able to quantify its diffraction efficiency as a function of wavelength and incident angle We will 
develop a general formula here for a uniform grating's diffraction efficiency, which we define as the fraction 
of transmitted radiant energy that is in the the first order. (Surface reflection losses are not included in this 
efficiency factor, but the reflection loss would be essentially equivalent to that of a conventional refractive 
lens surface and may be determined from the Fresnel reflection formulas. Standard means such as 
antireflection coatings may be used to mitigate reflection losses at the grating surface.) 

FIG. 6 illustrates the system of coordinates that will be used in the following analysis. The grating profile 
is defined by the relation y = f(x), where 
f(x) = x tana for |x| < | a/21; (18) 
f(x + A) = f(x). (19) 

(As discussed previously, the grating period a may be formally defined as a negative quantity ) The 
refractive index n has the form 



n = / n Q for y < £< x) 

\ n Q ' for y > f< x > . (2Q) 

For the present analysis, we assume that the incident beam is collimated with its propagation direction 
parallel to the x-y plane (i.e., perpendicular to the grating lines), and we assume that the beam's polarization 
vector is perpendicular to the x-y plane (transverse electric polarization case). We will represent the electric 
field in terms of its complex scalar amplitude distribution E, with the time dependence exp(-i(27r \)ct) being 
implicit, (t is the time coordinate, c is the free-space speed of light, and X is the free-space wavelength) 
The electric field in the space y < f(x) has the form 
E = A, exp(i<t>j) ( + backscattered component) for y < f(x) (21) 

where A, is the incident field amplitude (A s is constant) and *, is the incident field's phase distribtion which 
has the form 

$,| x>y = 2-n (x sinflo + y cose 0 ). (22) 

where do is the beam's incident angle at the grating substrate (FIG. 6). (The "backscattered component" in 
equation (21) represents reflected diffracted orders, which are not of interest for the present analysis ) 

Generally, a dispersion-compensating grating's period A would be at least one or two orders of 
magnitude greater than the wavelength X, so the electric field in the immediate vicinity of a grating facet can 
be represented fairly accurately by the form that it would take in the limit that A approaches infinity In this 
limit the transmitted field has the form of a collimated beam whose propagation angle 0, (FIG 6) is 
determined by applying Snell's Law at the grating surface: 
n Q ' sin(0, + a) = n 0 sin(0 o + a). (23) 

This representation of the transmitted field would not be valid in the immediate vicinity of the facet 
boundaries where the field's plane-wave form would be perturbed by edge scattering effects analogous to 
the shadowing and blocking effects that occur in the vicinity of a conventional Fresnel lens's facet sidewalls 
But these effects would generally be negligible due to the grating's extremely low-profile facet structure 
(For a dispersion-compensating grating, the facet angle a would typically be in the range of a couple 
degrees or less.) The transmitted field's plane-wave representation based on Snell's Law (23) would also 
not be valid in regions far from the grating surface where the aperture diffraction patterns from many grating 
facets would overlap, but we will only assume that the representation holds out to the level of the facet 
peaks at y = |A/2| tan*. We will represent the transmitted electric field in the immediate vicinity of a 
patncular grating facet approximately as 
E - A, exp(i* t ) for |x| < |a.2| and f(x) < y < |A2| tana (24) 
where A t the the transmitted field's amplitude and *, is its phase distribution, 

*i|x. y = 2it ~- (x sine, + y cose,). (25) 

(The reflected field in the immediate vicinity of a grating facet can similarly be represented approxi- 
mately as a collimated beam whose propagation direction is determined by applying the law of reflection at 
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the grating surface. The transmitted field amplitude A t and reflected amplitude can be determined by 
applying continuity conditions on the electric field and its normal derivative at the grating surface. For 
purposes of this analysis, however, we are not interested in analyzing reflection effects but are rather 
concerned with determining how the transmitted energy is allocated between diffracted orders, so we will 
leave A, undetermined.) 

As a consequence of the grating's periodic symmetry, the transmitted electric field outside of the 
grating region can be represented by a Rayleigh expansion of the form 
E = £ A m exp(i* m ') for y > |A/2| tana. (26) 

A m ' is the m-th transmitted orders amplitude and * m is its phase distribution. 



*m'|x.y = 2* J[ X- < x Slnd "' + V cos *™'>- < 27 > 

where 0 m ' is the m-th order propagation angle, which is determined by the grating equation, 
n 0 ' sina m ' = n 0 sine 0 + m\/A. (28) 

The total transmitted flux impinging on the plane y = |A/2| tana is proportional to |A t | 2 d4» t /dy, which \s 
is equal to |A t [ 2 (2tt n 0 '/x) cos0 t ; and the flux leaving this plane in the m-th order is proportional to |A m | 2 d$ m 
dy, or equivalents |A m '| 2 (2n n 0 Vx) cos0 m \ Hence the m-th order diffraction efficiency ij m (ratio of flux in the 
m-th order to the total transmitted flux) is 

20 n m = c °** m * |£^-| 2 . <29) 

The Rayleigh coefficients A m ' can be estimated by equating expressions (24) and (26) at y = |A'2| 
25 tana: 

A, exp(i*,) =» I A m exp(i* m ) at y = |A.2| tan*. (30) 

Multiplying both sides of this expression by exp(-i<f>o> ) and integrating from x = -A2 to x = A,2, we obtain 
an expression for A,' that can be substituted into (29) to obtain 

30 - coa fo* S1 nc2( * 2^1 A <sin*t - sm^' ) > . c31) 

The Rayleigh coefficients can alternatively be estimated by applying the continuity condition on dE'dy 
35 at y = |A/2| tana: 

|y (A t exp(i*,)) « §y = /A m ' exp(i* m ') at y = |A/2| tana. (32) 

Starting from this condition we obtain a diffraction efficiency estimate which is similar to (31) except that the 
cosine factors are interchanged: 

40 

f>m - l—l^ sinc2<tr A (sm^t - sin**')). (33) 

COS JS m K 

As a consequence of the simplifying approximation that we have made in representing the transmitted 
field by equation (24) in the grating region, the efficiency estimates (31) and (33) are not precisely 
consistent; but when |A| » x the sine 2 term in these expressions will only be of significant magnitude when 
& x is approximately equal to 0 m \ in which case the cosine terms will be practically equal. Hence, the 
difference between the cosine terms can be neglected and we can represent the diffraction efficiency as 

50 n Q ' , 

rj m = sinc 2 (7r ~-JT A (sine, - sine m )). (34) 

7, m is maximized and equal to 1 when the m-th order "blaze condition", $ m = e Xt is satisfied (i.e., when the 
m-th order diffraction angle 0 m ' defined by the grating equation (28) coincides with the angle of refraction 0, 
defined by applying Snell's Law (23) at the grating facet surface). 

By eliminating the sin0 m ' term between (34) and (28) we obtain the following alternative expression for 



n m = sinc 2 (*( m »« ( 35 > 
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where X a is the "blaze wavelength" defined by 
X B = A (n Q sinfl, - n 0 sm6 0 ) 36) 

(The sign of A would conventionally be defined so that X B is a positive quantity.) For a fixed diffraction 
geometry X B is constant and the efficiency in the first order is maximized when X = X B . It follows from (35) 
5 that the efficiency formula (34) satisfies the energy conservation relation, 
£ Dm = 1. (37) 

A dispersion-compensating gratign would not necessarily need to operate in the 1st (or -1st) diffracted 
order. According to equation (35), high diffraction efficiency in the m-th order can be obtained at any 
particular design wavelength X by choosing a blaze wavelength \ 8 = mX. But the diffraction bandwidth 
w would be very low if a high diffracted order is used, so a practical dispersion-compensating grating would 
need to operate in the 1st (or -1st) order. 



Conical diffraction geometries 

15 

The above analysis applies only to planar diffraction geometries in which the incident beam's 
propagation direction is perpendicular to the grating lines. (The diffraction efficiency derivation also 
assumes transverse electric polarization.) For the more general case, we will again consider the diffraction . 
20 geometry illustrated in FIG. 6, but we consider the incident beam's propagation direction to be inclined to 
the plane of the figure (i.e., the plane normal to the grating lines) by an angle 0. We introduce a third 
coordinate z with the z-axis running parallel to the grating lines, and we represent the electric field vector E 
in the incident medium as 

E = A, exp(i<fg ( + backscattered component) for y < f(x) (38) 
25 where A, is the incident field's amplitude vector and is its phase distribution, which has the form 
no 

^Ix.y.z = 2it x ((x sinflo + y cosflo) cos* + z sin*), (39) 

where do is the incident beam's projected propagation angle (i.e., the angle between the the substrate 
normal, or y-axis in Fig. 6, and the projection of the beam's propagation direction in the plane of FIG. 6). 
30 The transmitted field in the immediate vicinity of a particular facet has the approximate form of a 

refracted plane wave, 

E - A t exp(i<f> t ) for |x| < [A,2| and f(x) < y < \a.2\ tana (40) 

where A t is the transmitted field's amplitude vector and $ t is its phase distribution, 

35 *t|x. y ,2 = 2tt ((x sin0, + y cosfl,) cos*' + z sin*'). (41) 

The inclination angle * of the transmitted field's propagation direction relative to the plane of FIG. 6 and its 
projected propagation angle d x in the plane of the figure are determined by the following generalized Snell's 
Law relations: 

n Q sin* = n Q sin*, (42) 
40 n Q sin(0, + a) cos*' = n Q sin(0 o + a) cos*. (43) 

Outside of the grating region the field has the form of a generalized Rayieigh expansion, 
E = £ A^' exp(i<f> m ) for y > |A/2| tana, (44) 

where Arn is the m-th diffracted order's amplitude vector and <f> m ' is its phase distribution, 

45 Ix. y .z = 2tt -77- ((x sin0 m ' + y cos0 m ')cosa' + z sina). (45) 

The diffracted rays form a conical fan defined by the orders' common inclination angle *'. The m-th order's 
projected propagation angle * m is defined by the generalized grating equation. 
n Q sin0 m ' cos*' = n 0 sin0 o cos* + mX/A. (46) 

(Conical diffraction geometries would not typically be encountered with dispersion-compensating gratings 
50 since the dispersion of the diffracted orders' inclination angle *' is determined soley by Snell's Law (42) 
and hence cannot be influenced^ or corrected by the grating. But in some optical configurations the 
refraction-induced dispersion of *' at the grating surface might be offset by dispersion in another optical 
element, so a dispersion-compensating grating could operate in a conical diffraction mode.) 

An analysis similar to that outlined above yields the following general formula for the m-th order 
55 diffraction efficiency 

im = sinc 2 (7r - n ^> A (sins, - sin0 m ') cos*'); (47) 
or equivalently, 
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, m = sinc 2 (*( - m)), (48) 

where the blaze wavelength X B is defined by 

X B = A <n 0 ' sinfl, cos<*>' - no sine 0 cos<fr). (49) 

Nonuniform grating design procedures 



w The above analysis only considers diffraction of a collimated beam from a uniform grating structure. 
This model can fairly accurately represent a nonuniform grating's diffraction characteristics over any small, 
localized region on the grating; but for general design purposes it would be more convenient to use a 
design procedure that applies directly to nonuniform structures, rather than developing the design by 
piecing together a continuum of localized uniform grating geometries. We will develop such a generalized 

75 design procedure here, using an adaptation of the phase-matching principle. 

To introduce the procedure, we will first consider how the phase-matching principle can be applied to 
the simple problem of designing a conventional refracting lens surface. The incident and refracted beams' 
design geometries will be defined in terms of their respective phase distributions * and * , which are 
functions of positional coordinates. The beams* geometric wavefronts are surfaces of constant phase, and 

20 the geometric rays are trajectories that are tangent to the phase gradient (i.e.. normal to the wavefronts). 
Referring to FIG. 7, the incident beam's phase distribution $ defines, for any point P in the vicinity of the 
lens surface 22. the optical phase distance *|P (i.e.. 2-n x number of wavelengths) from a designated 
reference wavefront 23 in the incident beam to P. The refracted beam's phase distribution <* similarly 
defines the optical phase distance *'|P from a designated reference wavefront 24 in the refracted beam to 

25 P. The sign of a phase function will, by convention, be chosen so that the function increases in the beam's 
direction of propagation; hence will be negative in the vicinity of the lens surface. In terms of the 
geometric distances L and L (see FIG.7), the optical phase distances $|P and *'|P are defined as 

*|P = 2*-^-, , *'|P = -2*-^ (50) 

go where X Q and X 0 ' are the wavelengths on the lens's incident and exit sides, respectively. (Equations (2) 
define \ Q and x 0 ' in terms of the free-space wavelength X and the incident and exit optical media's 
respective refractive indexes n D and n G '.) Both phase distributions are extended across the lens surface 22 
to cover the spaces on both sides of the surface (as indicated by the dashed wavefronts in FIG. 7). 

The phase distributions are nonunique to the extent that changing either distribution by adding a 

35 constant merely shifts the position of the reference wavefront without affecting the structure of the beam 
represented by the distribution. If the optical surface in FIG. 7 is a component of an optical system, the 
phase distributions * and * may traverse other optical elements in the system, and the reference 
wavefronts may be located outside of the system. (In a region where the beam crosses itself, as in the 
space in front of a reflective surface, the phase distribution would be multiple-valued.) In general, the value 

40 of $ or <f>' at a point P would be calculated by summing the optical phase lengths, along a ray section 
joining the reference wavefront and P. of the ray segments in the different optical media traversed by the 
ray. (If the ray intercepts a grating element, the calculation must also include the discontinuous phase jump 
that occurs at the point where the ray intercepts the grating. A method for calculating the phase 
discontinuity across a transmission grating will be outlined later.) 

45 For any general point P located on the lens surface 22 in FIG. 7. the ray passing through P will traverse 
a total phase distance of (<t» - *>')|P between the reference wavefronts 23 and 24 (where the notation "(<f> - 
*')\P n indicates the value of the function <f> - *' at point P.). According to the phase-matching principle, this 
phase distance should be the same for all ray segments joining the wavefronts. so the function * - <t> 
should be constant over the entire lens surface: 

so (* - * )|P = C for any point P on the lens, (51) 

where C is a fixed constant Thus, the lens shape can be determined by specifying the location of any 
particular point P Q on the lens surface and stipulating that the phase function * - * at any other surface 
point P should be equal to its value at P Q ; i.e.. the surface is defined by the /elation 
(* - *')|P = (4> - 4>')|P 0 for any point P on the lens, (52) 

55 The phase-matching principle can also be applied to the design of a diffraction grating. Again, the 
incident and diffracted beams are specified in terms of their design phase distributions $ and * (for a 
particular design wavelength) and we initially suppose that the grating's substrate geometry is also 
specified. (The "grating substrate" is an imaginary smooth surface tangent to the grating grooves. The 
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surface is "smooth" in the sense that over any wavelength-scale region the surface is substantially flat.) As 
m the previously-discussed case of a uniform grating, the design is performed in two steps: First we 
partition the grating substrate into facet zones in such a way that the phase function $ - <J>' on each zone 
boundary is constant over the boundary and differs by 2tt between adjacent boundaries. (For a uniform 

5 grating, this condition determines the grating period A.) The facet zone geometry determines the diffracted 
beam's geometry, and the specified zone structure generates the specified beam geometry in the primary 
diffracted order (which would be a 1st or -1st order - an m-th or -m-th diffracted order of the specified 
geometry would be generated if the phase difference between adjacent zone boundaries were 2irrr\.) Having 
defined the facet zone structure we then define the facet surface geometry by specifying that each zone 

ro boundary defines a border of a facet surface over which * - #' is constant. (This condition defines the facet 
angle a for the uniform-grating case.) The facet surface shape determines the energy distribution between 
diffracted orders, and this specification maximizes efficiency in the primary order. 

The grating design is illustrated in FIG. 8. which shows a magnified portion of the grating surface 25. 
We will define the grating substrate surface 26 by a relation of the form 

75 f|P = 0 for any point P on the substrate. (53) 

where f is a specified function of positional coordinates and f|P is its value at point P. The first step in the 
design is to partition the substrate into facet zones (e.g., zone 27 in FIG. 8). Each facet face 28 will intercept 
the substrate over one of the zone boundaries 29. The facet zones are defined by the condition that the 
phase function <f> - $ is constant on each zone boundary and differs by 2tt between adjacent boundaries; 

20 hence the value of <t> - * on each zone boundary will be 2?rm + C where m is some integer associated 
with the boundary and C is a constant. The integer m may serve as an identifying index labeling the zone 
boundaries and associated intercepting facets, so the specification of the grating zone geometry can be 
stated 

f|P = 0 and (4> - <f>')|P = '27rm + C for any point P on boundary m. (54) 

25 The two relations in (54) define two surfaces (as illustrated by the dashed lines in FIG. 8). the intersection of 
which defines boundary m. (Any arbitrary value may be chosen for C in the above specification.) 

Each facet sidewall 30 (FIG. 8) would intercept a zone boundary 29 and may have a roughly straight 
profile transverse to the substrate. The detailed shape and orientation of the sidewalls may be be 
determined by practical considerations or may be optimized by a more refined optical analysis, but if the 

30 grating facet angle is small the sidewall geometry would have no significant influence on optical perfor- 
mance. The only general requirement that we will impose on the sidewall geometry is that its profile 30 
should be very short in relation to the facet surface profile 28 (FIG. 8). 

Each facet face 28 (FIG. 8) is bounded by two adjacent sidewalls, and the shape of facet m's surface 
between these bounds is defined by the phase-matching condition, 

35 (<f> - <f>')|P = 2?rm + C for any point P on facet face m. (55) 

(The facet profile shape defined by (55) would generally be slightly curved due to the wavefront curvature in 
the incident and diffracted beams, but in order to simplify lens tooling an approximately equivalent straight 
profile tangent to the ideal design profile could perhaps be used, at least over the grating's high-frequency 
regions, without significantly compromising optical performance.) 

40 The above design procedure defines a grating structure that diffracts the incident beam efficiently into 

an exit beam of the specified form for a particular design wavelength. (The grating zone structure defined 
by (54) determines the diffracted beam geometry, and the facet surface shape defined by (55) maximizes 
efficiency in the primary diffracted order.) If the grating's substrate shape is not predetermined, the extra 
degree of freedom would make it possible to design the zone structure to achieve achromatic imaging of at 

45 least two wavelengths. We will first consider how this can be accomplished for two very close wavelengths 
X g and X g + 5X, where 5X is an infintesimal quantity. The phase distributions $ and <f> will be defined as 
functions of wavelength over some spectral band containing X g . The incident phase function <f> may 
incorporate chromatic dispersion effects induced by elements preceding the grating, and the diffracted 
phase function <f> may be defined to compensate for dispersion induced by elements following the grating; 

so thus, the grating may compensate for dispersion induced by other optical elements as well as refraction- 
induced dispersion at the grating element itself. 

Applying condition (54) at wavelength X g , we obtain 
f|P = 0 and <$ - $')|P,Xg = 27rm + C|X g for any point P on boundary m. (56) 

where the "constant" term C may in the present context depend on wavelength. (The subscripts "P" and 
55 "Xg" in (56) indicate the position P and wavelength X g at which the functions are evaluated.) For wavelength 
X g + ax we also obtain the auxiliary condition 

(<t> - $')|P,\g + SX a 27rm + C|X g + 5X for any point P on boundary m. (57) 

Taking the difference of (57) and the second equation in (56), dividing by 5X, and taking the limit 5X --> 0. 
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( .M d&' > i 3 dC , 

a\ ax Ip, Xg dx! \g 

for any point P on boundary m. (58) 
Equation (58) will be consistent with (53) if we define the substrate's defining function f by 

i I = (M - Ml) | - dC| . <59) 

1 I P ax ax I P, x^ dxl x g 

75 The constant term dC/d\|\ g in (59) can be determined by specifying the location of any particular point P D 
on the substrate and stipulating condition (53) at this point. Equation (59) then translates to 

20 Ip ai: ax Ip.xg c ax ax Mp OP x g - 

With P 0 specified, the substrate geometry is defined by (53), the grating zone boundaries are defined by 

(56), and the facet shape in zone m may be defined by applying condition (55) at wavelength Xg in order to 
25 maximize diffraction efficiency at the design wavelength: 

(<f> - *')|P,X g = 277m + C|X g for any point P on facet face m. (61) 

(Any value may be chosen for the constant term C|X g in (56) and (61).) 

In optimizing the grating geometry for a wide spectra! band, improved dispersion-compensating 

performance may be achieved by using two widely-separated design wavelengths X r and x b .. Applying the 
30 second of conditions (54) on the zone boundaries, we obtain the conditions 

(<fr - $')|P.x r = 2irm + C|x r for any point P on boundary m (62) 

for wavelength X r , and 

($ - $ )|P,Xb = 27rm + C|X b for any point P on boundary m (63) 

for wavelength x bf where C|X f and C|X b are distinct constants associated with wavelengths X r and X b , 
35 respectively. Equations (62) and (63) will be consistent with (53) if we define the substrate's defining 
function f by 

f| P = ((* - *')|P iXr - (<& - *')|P.X b ) - (C|X r - C|x b ). (64) 

The constant term (C|X r - Cjx b ) in (64) can be determined by specifying the location of any particular point 
P 0 on the substrate and stipulating condition (53) at this point, f then translates to 

40 f|P s «* - *')|P.Xr - (* - *')|P.Xb )-((*- *)|Po.Xr - (* - *')|Po.Xt>) (65) 

and the substrate geometry is defined by (53). With the above definition of f equations (62) and (63) 
become redundant, so the grating zone boundaries can be defined by the conjunction of equation (53) with 
either of conditions (62) or (63) (using any value for the constant term C|X r or C|X b .) 

To determine the grating's residual dispersion, we need to be able to calculate the diffracted beam's 
45 phase distribution *' for any general wavelength x. given a specification of the grating structure and the 
incident phase distribution The primary diffracted order's beam geometry is defined by the phase- 
matching condition, 

(4> - 4>')|P.X = 27rm + C|X for any point P on boundary m. (66) 

where C is an arbitrary function of wavelength. Equation (66) implicitly defines on the facet boundaries. 

5 o (In order to maintain consistency between the phase distribution *' defined by (66) and its predefined 
values at wavelengths X r and X b , the values of C|X for X = X r and X = X b should be consistent with the 
values C|X r and C|X b in (62) and (63).) can be defined over the entire substrate surface by smoothly 
interpolating between the facet boundaries. If the grating structure is defined by the two-wavelength 
dispersion compensation procedure described above. *'|P.X may be implicitly defined on the grating 

55 substrate by eliminating the 2*01 term between (66) and either (62) or (63), and extending the resulting 
equivalence over the entire substrate. For example, using (62) we obtain 
(* - $')|P,X - C|X = (* - $')|P,x r - C|X r for any point P on the substrate, (67) 
where *'|P,x r is predetermined. 
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(In contrast to the case of a uniform grating, there is no precise mathematical criterion for identifying 
and distinguishing a 1st or -1st diffracted order in a general nonuniform grating's diffracted field; but over 
any small region of the substrate over which the local grating geometry is substantially uniform, the phase 
distribution <*> defined above coincides with the primary order's phase. The process of "smoothly 

5 interpolating" <t> between facet boundaries is imprecise only to the extent that the description of the local 
grating geometry as being "substantially" uniform is also imprecise. No ambiguity arises with this 
description if the facet profile dimensions are small in comparison to the grating's aperture dimension and 
the curvature radii of the substrate, facet boundaries, and incident wavefront.) 

Once 4> is determined on the grating substrate, its value can be readily obtained throughout the 

w diffracted beam. Referring to FIG. 9. the direction of a diffracted ray 31 (of any particular wavelength \) 
projecting from a point P on the substrate 32 is parallel to gradient vector V$' at P. ("V" is the gradient 
operator with respect to positional coordinates. The gradient vector V$' is indicated as 33 in FIG. 9.) Since 
$ is determined over the entire substrate, the tangential projection 34 of V<s>' on the substrate can be 
determined; and its normal projection 35 can then be determined from the condition that the magnitude k of 
75 V4>|P,X is equal to 2* n,X, where n is the value of the refractive index in the medium into which the ray 
transmits: 

|(V$'|P,\)| = 2* n/X. t (68) 

The gradient vector V4> at P is then obtained taking the vector sum of its tangential and normal projections. 
The direction of this vector determines the diffracted ray direction at P; and by ray tracing from P we can 

20 then determine the ray's complete trajectory through the optical system. The difference between the field's 
phase values at any two points along the ray is equal to the phase length of the ray segment joining the 
points (plus the phase jumps at any grating interfaces intercepted by the ray); so given the value of <*>' at P. 
we can readily determine its value at any other point on the ray. 

(The above procedure for calculating 4>' defines a general method for ray tracing through a grating. In 

25 analyzing beam transmission through a general grating element, the incident beam's phase distribution <f> 
can be extended past the grating by defining <f> = $' in the diffracted beam. From (66) we see that * has a 
jump discontinuity of -(27rm + C|X) across a point P on zone boundary m; and the phase discontinuity at 
substrate points between zone boundaries can be determined by interpolation. The phase discontinuity 
across the substrate surface at a general substrate point P determines a relationship between the phase 

30 gradient vectors 7$ on P's incident and exit sides, from which the relationship between the incident and 
diffracted rays' propagation directions can be determined.) 

The above procedure for determining only requires a specification of the grating's zone structure - 
the facet profile shape only determines the diffraction efficiency in the primary order and need not be 
specified to determine the diffracted field's phase distribution. Once the zone structure is determined, we 

35 can define the facet profile geometry by choosing an appropriate blaze wavelength x Bt defining the 
diffracted beam's phase distribution at this wavelength, and then defining the facet surface shape by 
applying the phase-matching condition (55): 

(<£ - <I>')|P,X B = 27rm + C|X B for any point P on facet face m, (69) 
where the constant term C|X B coincides with its value in (66) at X = X B . 

40 In applying the two-wavelength dispersion compensation procedure described above, the design 
wavelengths x r and X b would generally be iteratively adjusted to minimize the grating's residual dispersion 
over a specified spectral band, and the blaze wavelength x B would then be adjusted to maximize some 
measure of the grating's band-averaged diffraction efficiency. (The uniform-grating diffraction efficiency 
formula (35) may be used to calculate the diffraction efficiency over any small-scale region on the grating.) 

45 Since different parts of the grating may not exhibit optimum performance for exactly identical values of 

x r . x b , and X B , we may generalize the above design procedure somewhat by allowing the design 
wavelengths x n x b , and X B to be weakly-varying functions of aperture position, which may be independently 
optimized on different parts of the grating. In applying this generalization, the desired output phase 
distribution would be defined over the full wavelength range spanned by x r and \ b . However, the degrees of 

so freedom available with a single grating element do not allow the output phase function to be predetermined 
at more than two wavelengths, so we must draw a distinction between the ideal output phase distribution, 
ideal, and the actual phase distribution, In lieu of (65) the substrate's defining function f would in general 
be defined by 

f|P = ((* - *' ideal)|P.(X r |P) -(<*>-*' ideal)|P;(X b |P)) - ((* - ideal)|P 0 ,(X r |P 0 ) - (* - ideai)|P 0 ,(X b |P 0 )). 
55 (70) 

where the positional dependence of X r and X b is explicitly represented. We may define the facet zone 
structure by a generalization of (62), 

- ideal)|P,(X r |P) = 2^rm + C r for any point P on boundary m, (71) 

14 
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where C r is a fixed constant. (We could replace the constant term C|X, in (62) with a weakly-varying function 
of position of the form C|(X r |P) where C is a smooth function of wavelength, but the extra degree of freedom 
gained with this generalization is already provided through the wavelength-dependence of X r , so we simply 
substitute a constant C r for C|\ r .) Once the substrate geometry and zone structure are defined, the 

5 diffracted beam's actual phase distribution <*>' is implicitly determined by (66). (Whereas the identity 
between the functions in (62) and (66) required that the term C|X r in (62) and C|X in (66) be consistently 
defined, no such consistency is required between (71) and (66) since * and $ ideal are distinct.) The facet 
surface geometry may then be defined by a generalization of (69) 
($ - $')|P,(X 8 |P) = 27rm + C|(X B |P) for any point P on facet face m. (72) 

w (C may be an arbitrarily defined function of wavelength, but it should be defined consistently in (66) and 
(72).) 



General specification 



/5 



The design procedures described above can be applied to design a grating element that functions 
within an optical system to minimize the system's chromatic dispersion, provided that the system's design 
configuration (excluding the grating itself) is already specified. In practice, more sophisticated design 

20 methods would generally need to be employed since the optical system's design configuration would not 
normally be predetermined, but would rather be simultaneously optimized along with the grating design. 
The design principles outlined above may, however, serve to specify the invention and distinguish it from 
other grating applications. For this purpose, we consider a dispersion-compensating grating that forms a 
component of an optical system whose design configuration is specified, and we use the design specifica- 

25 tions for any optical elements preceding or following the grating in specifying the grating element's form 
and function. 

The optical system may, in the simplest configuration, consist only of the grating element itself. In 
general, the grating need not constitute the only dispersion-compensating element in the, system. For 
example, the system could achieve a high level of dispersion compensation by using two or more 

30 dispersion-compensating gratings in series, in which case the following specification would apply indepen- 
dently to each such element. 

The grating is formed as a relief pattern on the surface of a refractive optical element. Here, a 
"refractive optical element" is a transmissive element that, in conjunction with the grating, performs a 
primary optical function distinct from the grating surface's secondary dispersion compensating function. 

35 wherein the primary function alone could be accomplished by substituting a conventional refractive surface 
for the grating surface. To clarify this functional distinction, we can define the following structural distinction 
between the grating and the refractive element on which it is formed (see FIG. 10): Considering a grating 
surface 36 that separates two optical media with refractive indexes n c and n G . we can define the grating as 
a thin-film element comprising the region of space between a substrate surface (37) in the index-n 0 medium 

40 and a substrate surface 38 in the index-n Q ' medium. The optical element on the grating's index-n Q side can 
be regarded as a separate conventional refracting element whose surface is defined by the substrate 
contour 37. (The grating may be sandwiched between two such elements.) The refractive element can be 
thought of as being physically separated from the grating by a thin (i.e., zero-thickness) gap in which the 
refractive index is n Q \ A ray that impinges on the grating, say from the index-n 0 medium, will be deflected 

45 by two means, first by purely refractive means as it enters the index-n Q ' gap. and then by purely diffractive 
means as it traverses the grating into the index-n 0 ' exit medium. Thus, the grating and the refractive 
element would, in combination, perform their primary optical function through both refraction and diffraction. 
(The grating would typicaJly contribute only very marginally to their combined optical power, so the primary 
optical function alone could alternatively be achieved by eliminating the grating and slightly modifying the 

so refractive element's surface geometry.) In addition to performing this primary function, the grating would 
also perform a secondary dispersion compensating function (e.g.. balancing refraction-induced dispersion in 
the refractive element) via its own diffraction-induced dispersion. 

In the above context we define the grating surface's secondary "dispersion compensating function" 
generally as reducing chromatic dispersion in the optical system's output beam to a level significantly lower 

55 than what would be attainable if the grating surface were replaced by an optimally-configured conventional 
refractive surface. This definition distinguishes the grating's primary and secondary functions even when the 
refractive element on which it is formed comprises part of a conventional achromatic doublet whose primary 
function includes chromatic dispersion mitigation (e.g., the doublet may function to provide sufficient 
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degrees of freedom in the design to allow two-wavelength color correction, whereas the grating may 
function to provide sufficient additional degrees of freedom to achieve three-wavelength color correction) 

The optical system would function to channel a specified spatially coherent polychromatic incident 
beam into an output beam that is also substantially spatially coherent (to the extent possible within the 
s range of poss.ble design variations). The input beam is "spatially coherent" in the sense that any particular 
point ,n the port.cn of the beam's path preceding the system's first optical surface is traversed by only one 
(polychromatic) ray. Each such point will consequently also intercept a single incident wavefront and we 
will use one such wavefront as a reference in relation to which the input beam geometry will be defined 
In specifying the .input beam as being spatially coherent, we do not exclude the possibility that the 
w system may. under actual operating conditions, function to accept a spatially incoherent beam- but the 
system's dispersion-compensating mechanism will only be specified here in relation to a particular spatially 
coherent incident beam for which the system is intended to exhibit good dispersion-compansating 
performance. (For example, the dispersion properties of a system designed to focus direct sunlight may be 
characterized in relation to a perfectly collimated incident beam originating from a point at the center of the 
/5 sun disk.) 



A 
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polychromatic" beam comprises radiation of more than one wavelength w.thin a designated desiqn 
spectrum that may include as few as two wavelengths or may comprise a continuous band of wavelengths 
Oue to the effect of chromatic dispersion in the system, the beam's different chromatic components would 
generally diverge as the beam traverses the system, and the dispersion-compensating element functions to 
offset this effect in order to minimize chromatic dispersion in the output beam. In practice, the output beam 
will generally exhibit some residual dispersion, but we will describe the system in terms of a specified 
dispersionless (i.e.. spatially coherent) polychromatic output beam into which the system would ideally 
channel the input beam if chromatic dispersion could be eliminated. (Chromatic dispersion in the output 
beam would serve no intended purpose, and system's optical performance would be optimized if residual 
dispersion were absent in the output beam and the actual output beam conformed to the specified ideal 
beam geometry.) A particular wavefront following the system's last optical element will serve as a reference 
in relation to which the ideal output beam geometry will be defined. 

The input beam will be specified in terms of its phase distribution *. which is a function of positional 
coordinates and wavelength. For any particular wavelength X within the design spectrum, the set of points P 
satisfying the relation *|P.x = 0 defines a surface that coincides with the input reference wavefront By 
stipulating that this relation defines a common wavefront for all wavelengths X. we implicitly specify that the 
incident beam is spatially coherent. * will be defined throughout the region between the input reference 
wavefront and the grating surface. From a specification of the reference wavefront geometry * can be 
determined at any point within this region by ray tracing from the wavefront into the optical system The 
form of * in the v,cinity the grating depends on the design specification for all optical surfaces preceding 
the grating, and their chromatic dispersion characteristics are accounted for in the wavelength dependence 

The ideal output beam will be defined in terms. of its phase distribution *' ideal. The set of points P 
satisfying the relation * ideal|P.X = 0 (for any wavelength X within the design spectrum) defines the output 
reference wavefront. By stipulating that this relation defines a common wavefront for all wavelengths x we 
implicitly specry that the output beam is (ideally) dispersionless. *' ideal will be defined throughout' the 
region between the grating surface and the output reference wavefront. From a specification of the 
reference wavefront geometry, * ideal can be determined at any point within this region by reverse ray 
tracing from the wavefront back into the optical system. The form of *' ideal in the vicinity the grating 
depends on the configuration of all optical surfaces following the grating, and their chromatic dispersion 
characteristics are taken into account in the wavelength dependence of *' ideal. The output beam actually 
generated by the system will be defined in terms of its phase distribution *'. which is determined 
throughout the region between the grating surface and the output reference wavefront. *' can be calculated 
by ray tracing from the input beam through the grating. 

The grating structure is specified by first defining the substrate geometry, then defining how the 
subs rate is partitioned into facet zones, and then defining the facet surface geometry over each zone The 
substrate and zone geometries are defined so that the grating diffracts the incident beam into a primary (1st 
or -1st) diffracted order that substantially matches the specified form of the output beam and the facet 
surface geoemtry is determined to maximize diffraction efficiency in the primary order. 

The substrate geometry has a form that is determined so that the phase function * - *' ideal can be 
represented over the substrate substantially as the sum of some smooth function o that depends only on 
position, and some function C that depends only on wavelength: 
(* - * ideal)|P,X = 0|P + C|X 
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for any point P on the substrate and any wavelength X in the design spectrum. (73) 

(The sign in (73) and in other relations to follow means "substantially equal to", in a sense that will be 

discussed below.) Each zone boundary is defined as the locus of points P on the substrate over which n|P 

is equal to 2 7rm, where m is some integer that is associated with the boundary and may may serve as an 

identifying index labeling the zone boundary: 

Q|P = 2 for any point P on boundary m. (74) 

Given a specification of the substrate's shape and zone structure, the diffracted beam's phase 
distribution <f>' is implicitly defined on the zone boundaries by relation (66), 
($ . *')|P,x s 2 -rim + C|X for any point P on boundary m, (75) 

where the function C in (75) is chosen to correspond with the function C in (73). Interpolating between zone 
boundaries, *' is implicitly defined over the entire substrate by the relation 
(<t> - *')|P,X = Q\P + C|X for any point P on the substrate. (76) 

Taking the difference of relations (73) and (76), we see that the the specified grating structure generates a 
primary diffracted order whose phase error (*' - ideal) substantially vanishes over the entire substrate: 
($' - ideal)|P,X = 0 for any point P on the substrate and any wavelength X in the design spectrum: 
(77) 

and consequently the phase error substantially vanishes throughout the diffracted beam. 

If the grating substrate geometry could be chosen to satisfy condition (73) exactly, the diffracted 
beam's phase distribution would precisely match its ideal form ideal and chromatic dispersion would 
be completely eliminated in the optical system's output beam. There would not generally be sufficient 
degrees of freedom available in the choice of substrate shape to satisfy this condition exactly, however, so 
the specific sense in which condition (73) is to be interpreted requires some clarification. 

The substrate's design condition (73) can be equivalently stated by stipulating that the difference of the 
phase function ($ - ideal) between any two wavelengths Xi and X 2 in the design spectrum is substantially 
independent of position P on the substrate: 

(<£ - ideai)|P,Xi -(*-<$' ideal)|P,X 2 - C|Xi - C|X 2 for any point P on the substrate and any wavelengths \i 
and X2 in the design spectrum. (78) 

This follows directly from (73). It follows from (78) that the function (4> - <t> ideal)|P,X - C|X is substantially 
independent of wavelength X, and can hence be approximately equated to some function Q that depends 
only on position P over the substrate: 

(* - ideal)|P,X - C|X = Q|P for any point P on the substrate and any wavelength X in the design spectrum; 
(79) 

thus condition (73) also follows from (78). (The function C in (78) can be determined by arbitrarily defining 
its value C|Xi for any particular wavelength Xi, and then determining from relation (78) its value C|x 2 for any 
other general wavelength X 2 .) 

The substrate cannot generally be configured to satisfy relation (78) exactly, so the substrate shape 
would be chosen to minimize the positional sensitivity, over the substrate, of the expression on the left side 
of (78). One method by which expression (78)'s positional sensitivity can be minimized is to define the 
substrate by the condition that relation (78) should hold precisely for two particular wavelengths Xi and X 2 ; 
for example. Xi = X g and X 2 = X g + 5X, where X g is a chosen design wavelength and 5X is an infintesimal 
quantity. (One of the constants C|Xi or C|X 2 in (78) may be arbitrarily chosen, and the other may be chosen 
so that the substrate intercepts a specified point) From this design condition we obtain equations (53) and 
(59); and with the zone structure defined as in (56) the phase error (*' - <f>' ideal)|P,X g and its derivative with 
respect to wavelength X. (d($ - <f>' ideal)/dX)|P.X g , will vanish at any substrate point P. Hence, the diffracted 
beam's phase distribution *' will substantially match its ideal form <f>' ideal within a narrow spectral band 
including X g ; and in practice this design configuration would typically also exhibit good dispersion- 
compensating performance over a fairly broad spectral range. (If the optical system includes other 
dispersion-compensating elements, the substrate shape defined by (53) and (59) may also have the 
property that, at any substrate point P, the second and perhaps higher order derivatives, with respect to 
wavelengths X, of the function ($ - *')|P,X - C|X may vanish at X = X g ; consequently, with the zone structure 
defined by (56), the corresponding derivatives of the phase error (* - ideal)|P,X will also vanish at X = 
x g , resulting in further improved broadband performance.) 

For broadband applications, the substrate geometry may alternatively be defined by stipulating that 
relation (78) should hold for two widely separated design wavelengths X< = x r and x 2 = X b . From this 
design condition we obtain equations (53) and (64); and with the zone structure defined as in (62) the phase 
error ($' - ideal)|P,X will vanish at X = x r and at X = X b over the entire range of substrate points P. At 
intermediate and nearby wavelengths the phase error would generally be fairly minimal, resulting in good 
dispersion-compensating preformance over a broad spectral band. (If the optical system includes other 
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dispersion-compensating elements, the substrate configuration defined by (53) and (64) may also have the 
property that, at any substrate point P, the function (<f> - <t>' ideal)|P,X - C|\ is identical at three or more 
design wavelengths including X r and X b . With the zone structure defined by (62), the phase error (*' - *' 
ideal)|P.X would then vanish at all of these design wavelengths, resulting in further improved broadband 
5 performance.) 

Either of the above two design methods may be used to define the substrate geometry, or a more 
refined method may be used (e.g., we may define X r and X b as slowly-varying functions of positions over 
the substrate, as in (70)). From a general geometric perspective, these design approaches minimize 
expression (78)'s positional sensitivity over the substrate surface by defining the substrate geometry so that, 
w at any particular substrate point P, the substrate surface is approximately perpendicular to a space curve 
spanned by the gradient vector V($ - ideal)[P,X as X ranges over the design spectrum. (The gradient of 
the left-hand side of (78) is equal to the vector difference of two points on this space curve, and is hence 
approximately perpendicular to the substrate at point P; so the left-hand expression in (78) will exhibit 
minimal positional sensitivity over the substrate.) 
75 Any reasonable design approach used to minimize expression (78)'s positional sensitivity over the 

substrate will result in a grating design that dramatically reduces the optical system's chromatic dispersion 
in relation to a system that uses a conventional refractive lens surface in place of the grating surface, so in 
specifying the most general form for the grating substrate geometry, the precise design method employed 
need not be explicitly defined. Significant reduction in chromatic dispersion will be achieved if the chosen 
20 substrate geometry has the characteristic that, for any wavelengths Xi and X 2 in the design spectrum, the 
function ($ - <£ ideal)|P,Xt - (* - <f>' ideal)|P.X 2 in (78) is significantly less sensitive to position P over the 
substrate surface then it would be over an optimally-configured conventional refractive surface designed to 
replace the grating surface, so this general characteristic may serve to specify the range of allowable 
substrate geometries. (The function's degree of positional sensitivity over a surface can be quantified in 
25 terms of the magnitude of the tangential projection, onto the surface, of the function's gradient vector.) In 
stipulating that the function has sigificantly less positional sensitivity over the substrate than it would over a 
conventional refractive surface, it is to be understood that this comparison holds generally between any two 
points, one on the substrate and the other on the refractive surface, that are intercepted by a common 
incident ray (i.e.. the tangential projection of the function's gradient vector onto the substrate at a substrate 
30 point P would be substantially smaller in magnitude than the gradient vector's tangential projection onto the 
conventional refractive surface at a surface point that is intercepted by the same incident ray that intercepts 
P). The comparison would not hold in general between two points that are intercepted by different rays. (For 
example, residual dispersion at the peripheral zones of a dispersion-compensated lens would generally 
exceed the dispersion produced over the axial zone of an unconpensated lens; hence, the function's 
35 positional sensitivity over peripheral regions of the grating substrate would .typically exceed its sensitivity 
over axial regions of a conventional refractive surface.) 

With the substrate geometry specified so that the left side of relation (78) exhibits minimal or only 
moderate positional sensitivity over the substrate, the function C on the right can be defined by arbitrarily 
choosing its value C|Xi for any particular wavelength Xi and then defining its value C|X 2 for any other 
40 wavelength X 2 so that condition (78) holds substantially. (The error in (78) may in practice be minimized by 
equating CjXt - C|X 2 to the average or median value, over the substrate, of the left side of (78).) 

The grating's zone structure is defined by equation (74), where the defining function Q may be defined 
by stipulating relation (79) at some particular design wavelength X. For example, substituting equation (79) 
in (74) with the design wavelength X g substituted for X, we obtain the second equation in (56); and 
45 substituting the design wavelength x r for X. we similarly obtain (62). (The design wavelength x r may be a 
slowly-varying function of position over the substrate, as in (71).) Disregarding the specific design method 
used to define fl. it is evident from condition (78) that the left side of relation (79) will be substantially 
independent of position P over the substrate, and it will hence be possible to choose some wavelength- 
independent function Q that substantially satisfies condition (79). Any such function 0 may be used in (74) 
so to define the zone geometry. (The error in relation (79) originating from the wavelength dependence of the 
left side of the relation gives time to residual chromatic dispersion in the optical system's output beam, and 
this error may in practice be minimized by defining Q|P as the average or median value, over the design 
spectrum, of the left-hand expression in (79).) 

Each zone boundary borders an adjacent facet face, the face bordered by boundary m being 
55 designated as "face m". Each facet face extends over a substrate zone between two adjacent zone 
boundaries, with face m being bordered on one side by boundary m and on the other side by a short-profile 
sidewall joining face m to an adjacent facet face. (The sidewall's profile dimension is very small in relation 
to the facet face profile.) The surface shape of face m is substantially of a form defined by equation (72), 
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($ - $')|P,(X B |P) = 2 nm + C|(X B |P) for any point P on facet face m, (80) 

where X 9 is the "blaze wavelength which may in general be a weakly-varying function of aperture position. 
(The function is "weakly-varying" in the sense that it is substantially constant on any small region on the 
grating over which the grating zone structure is substantially uniform.) In practice the facet surface shape 
need not have the exact form defined by (80). For example, mold tooling operations may be simplified by 
using a straight-profile facet shape tangent to the siightly curved profile defined by (80); or the profile shape 
may deviate slightly from the form defined by (80) in order to achieve marginally improved diffraction 
efficiency over a broad spectral band. 

In a typical application, a dispersion-compensating grating would be formed on the surface of a Fresnel 
lens comprising large-scale Fresnel facets on which small-scale grating facets are superimposed (as 
illustrated in FIG. 2). Each Fresnel facet may be regarded as a separate lens element (e.g.. an annular lens 
in the case of an axially symmetric element), the configuration of which is defined by the above 
specification. 

A dispersion-compensating grating need not extend over the entire lens surface on which it is formed. 
Generally, a conventional refractive lens does not exhibit much chromatic dispersion near the lens axis, so a 
disperion-compensating grating would typically only be required over the lens's peripheral zones. 



Claims 

1. A dispersion-compensating transmission grating, formed as a relief pattern on the surface of a 
refractive optical element, that operated within an optical system to minimize chromatic dispersion, wherein 
the system comprises one or more optical surfaces (including the grating) that function conjunctively to 
channel a beam of polychromatic radiation comprising a distribution of wavelengths within a designated 
design spectrum, wherein 

the beam's phase in the portion of the beam preceding the grating may be represented by ..a function 
and its phase in the portion following the grating may be represented by a function 4> , wherein, 
$ and 4>' are functions of position and wavelength (the value of * or * at a general positional point P and 
wavelength X being indicated as "<f>\P,\ n or "4>|P,x". respectively). 

the relation <*|P,X = 0 defines, for any particular wavelength X, a set of points P comprising a wavefront 
surface that precedes the first surface intercepted by the beam and is independent of wavelength (i.e.. the 
incident beam at the system's first surface is spatially coherent), and 

in the system's intended mode of operation <J>' is intended to substantially match an ideal phase distribution 
4>' ideal, wherein the relation ideal|P,X = 0 defines, for any particular wavelength X. a set of points P 
comprising a wavefront surface that follows the last surface intercepted by the beam and is independent of 
wavelength (i.e., the exit beam at the system's last surface is intended to be substantially dispersionless). 
the grating substrate has a form specified by the condition that, for any two wavelengths Xi and X 2 in the 
design spectrum, the function 
($ - ideal)|P,Xi -(*-$' ideal)|P,X 2 

is significantly less sensitive to position P over the substrate surface than it would be over an optimally- 
configured conventional refractive surface designed to replace the grating surface, whereby the function 
satisfies an approximate relation of the form 
($ - $' ideal)|P,Xt ideal)|P,X 2 = C|Xi - C|X 2 , 

where C is some function that depends only on wavelength, and hence the term - * ideal)|P.X - C|X is 
substantially independent of wavelength X and can be represented substantially in terms of some function Q 
that depends only on position P: 
(* - ideal)|P.X - C|x = Q\P. 

the grating substrate is partitioned into a number of zones, each zone boundary being associated with a 
unique integer index m, wherein zone boundary m comprises the set of substrate points P for which the 
equation 

fl|P = 2 Trm , 
holds, whereby the grating generates a primary (1st or -1st) diffracted order whose phase function # is 
implicitly defined at any wavelength X and substrate point P by the equation 
(* - $')|P.X = Q|P + CjX, 

from which it follows that $' and ideal are substantially identical over the substrate and throughout the 
primary diffracted order: 
<f> « * ideal, 
and 
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for each index m, zone boundary m is uniquely associated with a grating facet face designated as "face 
m'\ wherein each facet face extends over a substrate zone between two adjacent zone boundaries, with 
face m being bordered on one side by zone boundary m and on the other side by a short-profile facet 
sidewall joining face m to an adjacent facet face, and wherein the surface of face m between these borders 
5 comprises a set of points P defined by the relation 
<* ■ *')|P.(X B |P) - 2 *m + C|(\ B |P). 

where x B is a blaze wavelength (which may be a weakly-varying function of position P). whereby the grating 
region in the vicinity of a general substrate point P will exhibit high diffraction efficiency at wavelengths near 
X 8 |P- 

w 2. A dispersion-compensating grating of claim 1. wherein the optical system of claim 1 is a Fresnel 

lens, the refractive optical element of claim 1 is a large-scale Fresnel facet on which the grating is formed 
as a molded pattern of small-scale grating facets, and the grating functions conjunctively with the lens and 
with other similar grating elements formed on other Fresnel facets to focus a polychromatic incident beam 
and to minimize chromtic aberration at the focal image. 

75 
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